AGGREGATION OF RAPIDLY VARYING RISKS AND ASYMPTOTIC 

INDEPENDENCE 



On 
O 
O 
(N 

G 

Q\ 

5? 

-I— > 



ABHIMANYU MITRA AND SIDNEY I. RESNICK 



Abstract. We study the tail behavior of the distribution of the sum of asymptotically independent 
risks whose marginal distributions belong to the maximal domain of attraction of the Gumbel 
distribution. We impose conditions on the distribution of the risks (X, Y) such that P(X + Y > 
x) ~ (const)P(X > x). With the further assumption of non-negativity of the risks, the result is 
extended to more than two risks. We note a sufficient condition for a distribution to belong to 
both the maximal domain of attraction of the Gumbel distribution and the subexponential class. 
We provide examples of distributions which satisfy our assumptions. The examples include cases 
where the marginal distributions of X and Y are subexponential and also cases where they are not. 
In addition, the asymptotic behavior of linear combinations of such risks with positive coefficients 
is explored leading to an approximate solution of an optimization problem which is applied to 
portfolio design. 
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1. Introduction 

Estimating the probability that a sum of risks X + Y exceeds a large threshold is important 
in finance and insurance, and hence m uch applied probabilit y research has been dedic a ted t o 



this goal. Recent results are found in Albrecher et alj 2009 1. Kluppelberg and Resnick 20081 



IWang and Tana 20091 . lAsmussen and Rojas-Nandavapal 20081 ]. lAlink et alj 2004 | . Embrechts and Puccett 
[20061 ] , IKo and Tana [20081 ] . Approximating this probability helps us evaluate risk measures for in- 
vestment portfolios as well as estimating credit risk. 

The problem is reasonably well understood when risks have regularly varying marginal distribu- 
tions but another important large class of risk distributions is the maximal domain of attraction of 
the Gumbel distribution, denoted MDA(A), where 

A(x) = exp{-e~ x }, x £ R, 

and MDA(A) is the class of distributions F for which there exist a n > 0, b n G R such that 

(1.1) lim nil — F(a n x + b n )) = lim nF(a n x + b n ) = e~ x , x£l 

n— >oo n—Kxy 

[Resnickl . Il987l . page 38]. It is als o well known th at the risks having distribution in MDA(A) are 
rapidly varying, i.e. — oo-varying [Resnick! . Il987l . page 53]. Within the class of risks (X,Y) with 
marginal distributions F,G £ MDA(A), results on aggregation of risks are known when X and Y 
are independent. However, actual risks are often not independent and a somewhat weaker concept 
called asymptotic independence, allows risks to be modeled as dependent and is more practical in 
many modeling situations. Risks X and Y in a maximal domain of attraction are asymptotically 
independent if for all x = (x\,X2), 

(1.2) lim H n (a^ Xl + b^,a^x 2 + b^) = G 1 (x 1 )G 2 (x 2 ) 
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where H is the joi nt distribution of X and Y and both G\ and Gi are non-degenerate extreme 
value distributions de Haan and Ferreiral . 120061 . page 229] . There are also results on aggregation 



of risks in the absence of asymptotic indepe ndence where the analogue of (II. 2p holds but with a 



limit distribution which is not a product; see iKluppelberg and Resnicld [20081 ] . 



This paper considers the case where the risks X, Y are asymptotically independent with marginal 
distributions F,G G MDA(A). We also allow one marginal tail to be lighter and the distribution 
with lighter tail does not necessarily belong to the maximal domain of attraction of the Gumbel 
distribution. 

Within the class of vectors (X, Y) satisfying asymptotic independence and marginal distributions 
F,G G MDA(A), two prominent but very distinct behaviors have been observed. 

(1) First, suppose (X, Y) are two iid risks with common distribution F which is sub exponential 
and F G MDA(A). Then X and Y are certainly asymptotically independent and 

(1.3) *n P * + Y> 4 =2. 

x->oo P{X > x) 

So one possible behavior is that the sum has a distribution which is tail equivalent to the 
distribution of a summand. 

(2) Very different tail behavior is exhibited in Theorem 2.10 of lAlbrecher et al.1 [20061 ] . who 



exhibit a distribution of (X,Y), with X and Y being asymptotically independent and 
identically distributed with common distribution F G MDA(A), but 

.. P(X + Y>x) 

hm — — — — = oo. 

x-+oo P(X > x) 

In Section 2, we give a set of conditions on the joint distribution of (X, Y), guaranteeing behavior 
of the first sort, namely, 

/ s , P(X + Y>x) 
1.4 lim -\- ± -^ = 1 + c, 

x^oo P(X > x) 

where c = lim^^oo P(Y > x)/P{X > x), the limit being assumed to exist. If c G (0, oo), our 
conditions imply that X, Y are asymptotically independent and each belongs to the maximal domain 
of attraction of the Gumbel. When X, Y are identically distributed, (jl.3p holds. Under the further 
assumption of non-negativity of risks, the result is extended for the case of more than two risks. In 
Section 3, we provide examples of distributions which satisfy our conditions. The examples include 
cases where the marginal distributions of X and Y are subexponential and also cases where they 
are not. We also show one example which does not satisfy our conditions but yet exhibits the tail 
equivalence between the distribution of the sum and that of the summand. Thus, our conditions are 
only sufficient. In Section 4, we summarize asymptotic behavior of finite linear combinations of risks 
with non-negative coefficients. In Section 5, we suggest approximate solutions for an optimization 
problem which is related to portfolio design. The paper closes with concluding remarks and a 
brief summary of numerical experiments which give a feel for whether asymptotic equivalence is a 
suitable numerical approximation for exceedance probabilities of aggregated risks. 

2. Asymptotic tail probability for aggregated risk 

2.1. Asymptotic tail probability for the sum of two random variables. We give conditions 
guaranteeing (jl.4|) . The constant c satisfies c = lmu^ry, P(Y > x)/P(X > x) G [0, oo). When 
c G (0, oo), X and Y are called tail-equivalent [Resnickl . Il971bl ] and then our conditions guarantee 



that both the marginal distributions F,G G MDA(A) and X and Y are asymptotically independent. 
When c = 0, our result extends to the case where G, the marginal distribution of Y, does not belong 
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to the maximal domain of attraction of the Gumbel distribution and where X and Y need not be 
asymptotically independent. 

2.1.1. Assumptions. Suppose (X, Y) is a pair of random variables satisfying the following set of 
assumptions. 

(1) The random variable X has a distribution F whose right endpoint xq is infinite; that is, 

(2.1) xq = sup{x : F(x) < 1} = oo. 

Further F £ MDA(A) so that (11.11) is satisfied wi t h centering con stants b n £M and scaling 
constants a n > 0. Equivalently ( de Haanl [1970f |. iResnickl 19871 . page 28, 40-43] ) there 
exists a self-neglecting auxiliary function /(•) with its derivative converging to 0, such that 

(2.2) hm gft+jgW) = e -.. 

V ' *->oo F(t) 

(2) The random variables X and Y have distribution functions F and G such that 

hm ^ = ce [0,oo). 

z-+oo £ (x) 

(3) The conditional distribution of Y given X > x, satisfies for all £ > 0, 

lim P(|Y| >i/(x)|X>x) =0, 

x— >oo 

where /(x) is the auxiliary function corresponding to the distribution of X given in (12. 2ft . 

(4) and symmetrically assume for all t > 0, 

lim P(|X| >t/(x)|y >x) =0. 

(5) For some L > 0, suppose 

,. P(Y>L/(x),X>L/(x)) 

nm 7T7^ \ = °- 

2.1.2. The main result. The assumptions allow us to conclude aggregated risks are essentially tail 
equivalent to individual risks. 

Theorem 2.1. Under Assumptions^^ in Section \2.1.1l we have 

P(X + Y >x)~ (l + c)P(X >x), x^oo. 

2.1.3. Comments on the assumptions. Before giving a proof of Theorem 12. 11 we discuss implications 
of the assumptions. 

Remark 2.2. (1) When F g MDA (A), we ma y choose a n , b n appearing in (11.1 p as b n = bp{n), 
in = f(b n ). See (Resnickl . 1 19871 . page 40] or Ide Haan and Ferreiral 20061 ] . 
(2) If c £ (0, oo), then our assumptions guarantee both marginal distributions F,G £ MDA(A) 
and also that (X, Y) are asymptotically in dependent. From Assumption \T\ F £ MDA{A) 
and since F and G are tail-equivalent, from IResnickl 1971bl ] we get that G S MDA(X). For 
asymptotic independence, define, 

(2.3) b F (t) = inf{ S : —^00 > «> = (i^)"^)' 

and similarly bc{i). From [de Haan and Ferreiral . 120061 . page 229], if F, G & MDA(A) and 

(2.4) lim P(X > Mt),Y >bc(t)) _ 

v ! i-oo p(x > b F {t)) 
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then (X, Y) are asymptotically independent according to (|1.2p . When c £ (0,oo), Assump- 
tion [3] implies (|2.4p . To verify this, note first that Assumption [3] implies 

(2 . 5) lim PV>x.Y>x) < , im P(X> m ,Y>z) _ 

V ' x^oc P(X > X) - z^oo P(X > x) 

since f(x)/x — ► as x — > oo [Resnickl . 1 19871 . page 40]. If c > 1, then for sufficiently large t, 
b F (t) < ba(t) and therefore, using (|2.5|) . 

ppr>M*).^>M*)) < p(A>6 F (t),y>6 F (Q) 

t-™ P(X > b F (t)) ~ ti™ P(X > b F (t)) 

as required. A similar verification can be constructed for the case < c < 1. For c = 
l,b F (t) ~ ba(t). Hence, 

f(b F (t)) f(b F (t)) 



b G (t) b F (t) 

So, 



0. 



p(x>6 F (t),r>6 G (t)) < Hm p(x>fe F (t),r>/(M*))) 

A™ P(A > b F (t)) ~ t^™ P(X > b F (t)) 

= (by Assumption [3] and ((2H]) ) . 

(3) The auxiliary function /(•) can be replaced by any asymptotically equivalent function /(■); 
that is, if hnxj^oo f(x)/f(x) = 1, and if Assumptions El [H [5] hold with /(•), they also hold 
with /(•) replacing /(•) and vice versa. Since the mean excess function 

e(x) = E(X -x\X > x) 

is asy mptotically eq uivalent to any auxiliary function f(x) ( [Embrechts et all 1 199 71 . page 
143], [Resnicld . Il987l . page 48]), e(x) can also be taken as an auxiliary function. 

(4) If c = lirn^oo G{x)/F{x) = 0, we do not need Assumption |4] to conclude our result. 

(5) An easier proof of the result can be given if Assumption [5] holds for all L > 0. But here we 
provide an example to show the importance of the weak version of Assumption [5j 

Example 2.3. 

X = -\og{U), Y = -log(l-U), C/~ Uniform (0,1) 

It is obvious that in this case both X and Y have distribution Exponential(l). So, in this 
case, the auxiliary function is f(x) = 1. Choose L such that exp(— L) = |, and 

P{X > Lf{x), Y > Lf(x)) _ P{U < exp(-L), 1 - U < exp(-L)) 



P{X > x) P{X > x) 

A> _ 



P{\<U < 3 



OO. 



P(X > x) 2P(X > x) 

Therefore, this particular choice of L does not satisfy Assumption [5j The distribution of 
(X, Y) is a special case of Example E3] which discusses certain L which do satisfy assumption 

El 

(6) If, however, both X and Y are non-negative risks, and Assumption E] is strengthened to 
hold for all L > 0, then Assumptions El and U] will be automatically satisfied. The proof of 
this follows from lirn^oo f(x)/x = 0. 
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(2.6) 



(7) Similar limit r esults are found in Lemma 2.7 of lAlbrecher et al.l 20061 ] and Theorem 2.1 of 
IKo and Tana 20081 ] . They have assumed that one of the marginal distributions of the two 
asymptotically independent variables X and Y, say the distribution of X, is sub exponential, 
(i.e. X £ S, where S is the set of all sub exponential distributions) and worked on finding 
conditions for the tail-equivalence of the marginal distribution of X and the sum X + Y. 
Our assumptions are different: We assume that one of the marginal distributions of the 
two asymptotically independent variables X and Y, say the distribution of X, belongs to 
the domain of attraction of Gumbel, i.e. X € MDA(A). We do not assume the marginal 
distribution of X is sub exponential. 

In examples where the marginal distributions of the two asymptotically independent and 
identically distributed random variables X and Y belong to the class MD A(A) n S, an 
issue is the relative strength of our conditions versus those of Theorem 2.1 of lKo and Tang 
[20081 ] . We can not show either set of conditions implies the other. Below we present an 
example which sat isfies our set o f cond itions, but does not satisfy the set of conditions given 
in Theorem 2.1 of IKo and Tang! [20081 ] . Thus our set of conditions is not stronger. 



Example 2.4. Suppose, X = exp(-Xi), Y = exp(X2), where (X\,X2) is bivariate normal 
with correlation p € (0, 1). For simplicity, assume each Xi has mean and variance 1. It 
is well known that lognormal distribution belongs to the class MDA(A) n <S. In Example 
3.61 we show (X, Y) satis fy our set of conditio ns. Here we show that this example does not 
satisfy Assumption 2.1 of IKo and Tang) [20081 ]. i.e. for all x* > 0, 

P{Y >x-t\X = t) 



lim sup sup 

x— »oo x*<t<X 



oo. 



P(Y > x-t) 

From the exchangeability of X and Y, it is obvious that (|2.6|) holds even if the role of X 
and Y is interchanged. 



sup 

x*<t<x 



(2.7) 



P(Y > x - t\X = t) 
P{Y >x-t) 



sup 

x*<t<x 



= / log(a:-t)-plogt 

$ (log(s - t)) 



> 



ft I log(x/2)-p\o g (x/2) 

*(log(x/2)) 



<I> 



V^- 



log(x/2) 



oo. 



$ (log(x/2)) 

The inequality above follows from choosing x enough large so that x/2 > x* and putting 
t = x/2. The last convergence follows from the fact that nor mal distribution belongs to 
the class MDA(A) and hence $ is — oo-varying JResnicki . 1 198 71 . page 53]. Note, < p < 1 
entails 



V^- 



< 1. Hence, from (|2.7[) it is obvious that (|2.6|) holds. 



2.1.4. Proof of Theorem \2.1\ We prove Theorem 12.11 using a Proposition and a Lemma, which 
we prove first. Note, we do not need the assumption that the marginal distributions are sub- 
exponential, which is a necessary condition in the case where X and Y are independent. 

Proposition 2.5. Under Assumptions^ and [31 of Section \%.1.1\ we have 



(2.8 



lim P(Y < a n z\X > a n x + b n ) = 1 



{z>0}> 



and from Assumptions^ and\Q of Section \2.1.1\ we have 



(2.9) 



lim P(X < a n z\Y > a n x + b n ) = 1 



{2>0}> 



z^O, x G 



z / 0, i £ 
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Proof. The self-neglecting property of the auxiliary function /, i.e. 

r /(* + */(*)) i clB 

, lim 777^ = !> x€R, 

implies that 

(2.10) lim P(\Y\ > zf(t)\X > t) = lim P(\Y\ > zf{t)\X > t + xf(t)). 



t— »oo 



Hence, by noting that a n = f(b n ) and lim. 
part is proved similarly. 



t— *oo 

x> b n = oo, the result follows from (|2.10p . The second 

D 



Lemma 2.6. (i) Assumptions^^ and f51 of Section \2.1.1\ imply that the sequence of measures 

nP[a-\X-b n ,Y)£(dx,dy)] 

converges vaguely on ([— M, oo] x [—00,00]) as n ^ 00, to £fte /imif measure mi )00 (dx,dy) = 
e~ x dxeo(dy), for some M > L(from Assumption^ of Section \2.1.1\) such that —M is a continuity 
point of ~ n for all n. 
(ii) Assumptions^ 0, anrfrj] of Section \2.1.1\ imply that the sequence of measures 

nP[a-\Y-b n ,X)G (dx,dy)] 

converges vaguely on ([— M, 00] x [—00,00]) as n ^ 00, to £/ie /imii measure m2 )0 o{dx,dy) = 
ce~ x dxeo(dy), for some M > L(from Assumption^ of Section \2.1.1\) such that —M is a continuity 
point of - " for all n. 

Remark 2.7. Since all the discontinuity points of — — - for all n is countable, choice of such an 
M > L is not a problem. Moreover, the M in the two parts of the lemma (i) and (ii) may be chosen 
to be the same. 

Proof. We consider convergence of the measures evaluated on certain relatively compact regions 
which guarantee vague convergence. 

Region 1: (x, 00] x [— 00, y], x > -M,y / 0. As n —> 00, 

~X - b n 



nP 



X-b n Y 
> x, — < y 

(Xn (Xn 



nP 



> x 



d„ 



P 



Y < X-bn 

— < y\ > x 

(Xn. (Xn, 



-> e Xl {y>o} = mi )00 ((x,oo] x [-00, y]) 

by Proposition 12.51 

Region 2: [-M, x] X (y, 00], x > —M, y / 0. Since — M is a continuity point of x ~ bn for all n, 
as n — > 00, 



nP 



%r X-b n Y 

-M < < x,— > y 



nP 
nP 
nP 



, , X - b n Y 

-M < < x, — > y 

a n a n 

r^L^>_ M ,— >y 



X -b n Y 
> x, — > y 

(Xn. (Xn, 



X-br, 



> -M 



P 



nP 
> (e M 



X-b n 



> X 



(In 



— > y\ > -M 

a n a n 

Jr — > y\ > x 

(Xn Q"n. 



') l {y<o} = mi,ao([-M,x] x (y,oo\), 



by Assumption [T] and Proposition 12.51 
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Arguments for convergence on the following regions follow in a similar fashion using Proposition 



Region 3: (x, oo] x (y, oo], x > -M, y / 0, 

Region 4: [— M, x] x [-oo, y], x > -M,y^ 0. 
This concludes the proof of vague convergence on part (i) . 

The proof of part (ii) is similar, only notice that if c = 0, we do not need the Assumption [U In 
this case, note that the limit measure rri2,oo(dx, dy) is a zero measure. Also note, using Assumptions 
[H and El we get 

~Y-b n 



nP 



> -M 



ce M = 0, 



which is enough to prove the convergence in this case. □ 

This leads to a formal statement of the main result. 

Theorem 2.8. Under the Assumptions in Section \2.1.1[ 

P(X + Y > x) , 
2.11) lim -\- ± —L = (l + c). 

x^oo P(X > x) 

Proof. Choose M as in Remark [2771 We split P(X + Y >b n ) as 

P(X + Y >b n ) = P(X + Y>b n ,X >b n - Ma n ) + P(X + Y>b n ,Y>b n - Ma n ) 
- P(X + Y >b n ,X >b n - Ma n , Y > b n - Ma n ) 

(2.12) + P(X + Y >b n ,X <b n - Ma n , Y <b n - Ma n ). 
Using Assumption [1] and (|2.5p , we get 

nP(X + Y > b n ,X > b n - Ma n , Y >b n - Ma n ) < nP(X > b n - Ma n , Y > b n - Ma n ) 

(2.13) = nP(X > b n - Ma n ) p{J{ > - _ ^ >e .0 = 0, 

since b n — Ma n — > oo. Now, consider the convergence of the last term of (12. 12ft mutiplied by n. 

P(X + Y >b n ,X <b n - Ma n , Y <b n - Ma n ) < nP(X > Ma n , Y > Ma n ) 
(2U) P{X > Mf(b n ),Y> Mf(b n )) ^ P(X > Lf(b n ),Y > Lf(b n )) 

1 ' ' ~ P(X > b n ) ~ P(X > b n ) ~" ' 

by ()2.ip and Assumption [5J 

To deal with the first term of (|2.12p mutiplied by n, we first define a function T as T : [— M, oo] x 

[— oo, oo] h^ (oo, oo] by 

T ( x \ = } x + y > iiy> ~°°' 

[0, ify = -oo, 

and hence 

(2.15) nP{X + Y >b n ,X>b n - Ma n ) = nP{a~ l {X - b n , Y) e T^((0, oo]) n [(-M, oo] x {0}]). 

Note, that every set in the space [— M, oo] x [—00,00] is relatively compact, and hence so is 
r*~((0, 00]) n [(— M, 00] x {0}] = S (say). Also, since the limit measure mi )0O is concentrated on 
[-M,oo]x{0}, 

(2.16) mi i0O (<55) = m ltO0 (SS D [[-M, 00] x {0}]) = m ljOO ({0} x {0}) = 0. 
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([27L5]) and (1236]) . we get 
nP(X + Y>b n ,X>b n - Ma n ) -» m lt00 (S) 



Hence, using Lemma 
(2.17) 
Similarly, 

(2.18) nP(X + Y>b n ,Y>b n -Ma n )^m 2 ,oc 

Hence, using (|215]l . (|2T3l . (gUP , (I2T7D and IjZlBl) . we get, 

P(X + Y > x 



1. 



(5) 



lim 
and we conclude our result. 



P(X > x) 



lim nP(X + Y > b n ) = 1 + c, 



□ 



One immediate application of Theorem 12.81 is to t he subexponential fami ly of distributions de- 
noted S. The class MDA(A) f\S has been studied in [Embrechts et all 1 19971. page 149] and several 



sufficient conditions for belonging to this class are given in iGoldie and Resnicld [19881 ] . Corollary 
12.91 gives an additional sufficient condition and follows directly from Theorem 12.81 Example 13.21 
exhibits a distribution which satisfies the conditions of this Corollary. 

Corollary 2.9. Suppose, F £ MDA{K) with auxiliary function f(x) as described in Assumption 
d of Section \2.1.1\ Suppose, also, linx^oo f{x) = oo ; and for some L > 0, 

(2.19) !-P*tf = 0. 



lim 

x— >oo 



[F{Lf{x 
F{x) 



Then, for X and Y iid with common distribution F we have, as x — ► oo, 

P[X + Y >x] ~2P[X >x], 
and therefore, if F concentrates on [0, oo), F £ S. 

Following Remark 12 .2I |3|). it is enough to check (J2.19P with any f{x) satisfying f(x) ~ f(x). Note 
also it is n atural to add the as s umpt ion f(x) — ► oo, since if F G MDA(A) n S, then necessarily 



x 



Goldie and Resnicki . Il988l |. 



2.2. Asymptotic tail probability for the sum of more than two non-negative random 

variables. Suppose, among the risks Xi,X 2 , ■ ■ ■ Xd, there is no heavier tail than X\ in the sense 

that it is not true that 

Fi(x) 



lim 

x^oo Fi(X) 



oo, 



2,...,d. 



Assume X\ satisfies Assumption [T] of Section 12.1.11 and that X\ , X 2 , . . . Xj pairwise satisfy the 
Assumptions [3] and U] of Section [2.1.11 with the auxiliary function /(■) of X\. By this, we mean for 
all pairs 1 < i ^ j < d, and for t > 0, 

P(X j >tf(x),X i >x) 



lim 

x— >oo 



P(X t > X) 



0, 



P(X j >tf(x),X i >x) 



which implies 

(2.20) lim - V " J - ; J ^ / '" t 1 ~> = o. 

Also, suppose, the risks X\,X 2 , ■ ■ ■ X^ pairwise satisfy Assumption [5] of Section [2.1.11 with auxiliary 
function /(•) of X\ so that for 1 < i < j < d, there exists some L« > 0, such that either 

P(X i >L ij f(x),X j >L ij f(x)) 



lim 

x— >oo 



P{X t > X) 



0. 
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or. 

Hm P(X i >L ij f(x),X j >L ij f(x)) = Q 

x->oo P{Xj > X) 

In either case, we have, for 1 < i < j < d, for some L%j > 0, 

(2.21) lim — y — JJ ; h \ 3Jy n = 0. 

x^oo P{X\ > X) 

Under the additional assumption of non- negativity, Theorem 12.81 can be extended to more than two 
risks. 

Corollary 2.10. Assume, X±,X2, ■ ■ ■ Xj are non-negative random variables which pairwise sat- 
isfy Assumptions 3|j)[|o/ Section \2.1.1\ with the auxiliary function /(•) of X\. Moreover, the 
distribution of X\ satisfies Assumption^ of Section \2.1.1\ and suppose 

P(Xi > x) 
2.22 lim -±-±-—L = Ci e 0, oo , i = 2, 3, . . . , d. 

x^co F(Xi > X) 

Define, Sj = X\ + Xi + . . . Xj , 1 < j < d and we have, for x £ R, 

d 
(2.23) lim nP{S d > a n x + b n ) = (1 + V Ci)e~ x 

n — >nn * — ^ 



i=2 



and hence 



(2.24) lim P \ Sd > X \ = (1 + Y a) 

1=2 

Remark 2.11. (1) Asymptotic independence of the random variables: Suppose, for all i, Cj S 

(0, oo). Then for any 1 < i ^ j < d, the pair (Xj,X,) is asmptotically independent by 
Remark 12. 2l f2"jh Since the random va riables are pai rwise asymptotically independent, they 
are also asymptotically independent Resnickl . 1 198 71 . page 291]. 

(2) Non-negativity of random variables: The only additional assumption added to the list in 
Section 12.1.11 is that the random variables are non-negative. 

(3) Relaxation: We have shown in (|2.20p and (|2.2ip that pairwise satisfaction of Assumptions 
H H E of Section I2XT1 implies that for 1 < i ^ j < d, for t > 0, 

x->oo P(Xi > x) 

and for 1 < i < j < d, there exists Lij > 0, 

Hm P(X j >L ij f(x),X i >L ij f(x)) Q 

a-^cx) P(X\ > X) 

We will show that actually these conditions are enough to get the desired conclusion. 

Proof. We prove the result by induction under the relaxation Remark I2.11l |3|) . The base case of the 
induction for d = 2 is already proved in Theorem 12.81 so suppose, the result is true for d = k > 2 
and we have 

k 

(2.25) lim nP(S k > a n x + b n ) = (1 + V a)e~ x 

n— »oo ^— ' 

8=2 
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and 

(2-26) lim ^ > X \ = l+Ya 

V ' z^oo P(X 1 > x) ^ 



Therefore, we have 



P(X k+1 > x) c k+1 



(2.27) lim ;." +1 ^ = -^ £[0,oo). 



x->oo P(S k > x) i + sy 



:2 °« 



We will use Theorem 12.81 with X = S k and Y = X k+ \. It remains to check the Assumptions 
in Theorem 12.81 For Assumption [H note that S k is tail-equivalent to X\ and use the fact that 
F £ MDA(A) is closed under tail-equivalence. Assumption [2] is already checked in (J2.27J) . 

Note that from the induction hypothesis P(S k > x) ~ P(uf =1 (A"i > x)) and from the positivity 
of the risks [S k > x] D U^ =1 [Xj > x\. From these two facts, it easily follows that 

(2.28) lim P((5 t >,)nfe(X, >,))') = 

rr^oc P(5 fc > x) 

Since Sfe and X± are tail equivalent, by iResnickl [1971bl |. the auxiliary function /(•) of S k is 
asymptotically equal to the auxiliary function /(•) of X-y. Therefore, given e £ (0, 1), there exists T 
such that for all x > T, f(x) > ef(x). We now check Assumption [3l For t > 0, x > T, using (J2.28D . 
as x -> oo, 

P(\X k+1 \ > tf(x)\S k >x)< P(X k+1 > tef(x)\S k > x) 

_ P(X k+1 > tef(x),S k > x) P(X k+1 > tef(x),S k > x,uf =1 {X 4 > x}) 



P(S k > x) P(S k > x) 

tef(x),U* =1 {Xi > x}) < Zi=iP(Xk+i > tej 
P(Sk>x) ~ (1 + YH =2 *)P(X 1 >x) 



^ P(X k+1 > tef(x),ut =1 {Xj > x}) < Ei=i p ( x k+i > tef(x),Xj > x) Q 



by (|2T20l) . 

For Assumption HI if c^+i = 0, following Remark l2.2H pij). there is no need to check assumption [H 
So, suppose, c k+ i > 0. Then for any t > 0, as x — > oo, 

P(\S k \ > tf(x)\X k+l >x)< P(S k > tef(x)\X k+1 > x) 

k 

<J2P(Xi>tef(x)/k\X k+1 >x) 



k 



i=l 



P(X t > tef(x)/k,X k+1 >x) P(X 1 >x) 

P(Xy > X) P{X k+l > X) 



E r\y\j ^ it J {X)/K,y\ k+ l ^ X) J-^l ^ X) 
P(T, -» >A P( Y, , , -> -r\ ~^ ' 



For Assumption [U we know from the assumptions in the statement of Corollary 12.101 that the 
random variables satisfy Assumption [5] of Section 12.1.11 pairwise with auxiliary function /(•) of 
X\. Thus, for 1 < i < j < d, (|2.2ip holds. We check Assumption [5] with L = kL max /e, where 
Lmax = maxi<j<fc Lj 5 fc + i (recall, equation (|2.2ip ). Then, for sufficiently large x, using /(•) as the 
auxiliary function of S kl 

P(X k+1 > Lf(x),S k > Lf(x)) < P(X k+1 > Lef(x),S k > Lefjx)) 
P{S k >x) ~ P{S k > x) 
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< P{X k+ i > kL max f(x),U^ =1 {Xj > L max f(x)}) 

P(S k > x) 

< P(X k+1 > L iM1 f(x),u! =1 {Xi > L iM1 f(x)}) 

P(S k > x) 

< T,i=i p ( x k+i > L hk+1 f(x),Xj > L itk+ if(x)) 

P(S k > x) 

Eti HXk+1 > hh+ll^h x * > hh+llM , o 

by (|2.2ip . This completes the induction proof. D 

3. Examples 

This section shows a few of the many models that satisfy the Assumptions in Section 12.1.11 In 
all examples, both X and Y are non-negative random variables and it is straightforward to extend 
these examples to the d-dimensional case and show the assumptions of Corollary 12. 101 are satisfied. 

Our conditions are only sufficient and we exhibit one example where our conditions do not hold, 
but tail equivalence in Theorem 12.81 holds true. Finding a necessary and sufficient condition for the 
conclusion of Theorem 12.81 is still an open but subtle and difficult issue. 

Example 3.1. Suppose X\,X 2 ,X 3 are iid with common distribution F, where for a > 1, 

exp{— (logx) a }, if x > 1, 



F(a;)- 

1, it x < 1. 

Define, 

X = X 1 f\X 2 , Y = X 2 AX 3 

It is easy to check X and Y are identically distributed with the common distribution F\, where 

Fi(x) = exp(-2(logx) a ), x > 1. 
It can be checked that F\ is a Von-Mises function; that is, it satisfies, 



F X F' X 



1. 



m=^m= - „ * ! , *>i 



( F [r 

a sufficient condition for F\ G MDA(A), and 

Fi(x) _ x 

F[{x) 2a(\ogx) 

serves as an auxiliary function [Resnicki . 119871 . page 40]. Also, (|2.1|) is obvious and therefore, 
Assumption [1] of Section [2.1.11 is satisfied. Checking Assumption [2] is straightforward, so consider 
AsumptionEl Fix t > 0, recall f(x)/x — > 0, and note as x — ► oo, 

P(X >x,Y >tf{xj) P(X 1 >x,X 2 >x\/tf(x),X 3 >tf(x)) 
P(X > x) ~ Pjxl >x,X 2 > x) 

P(X 1 >x,X 2 >x,X 3 >tf(x)) r, fY ^. ff ^ n 

P(X 1 >x,X 2 >x) = P{Xi > tf(x)) "> °' 

since f(x) — > oo. Assumption S] is verified the same way. For Assumption \5\ we have with L = 1, 
P(X>/(s),y>/(s)) P(A 1 >/(x),X 2 >/(x),X 3 >/(x)) 
P(A>x) P(X X >x,A 2 >s) 



12 
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HfWY 



F{ x y 

exp < — 2(logx) a 



exp{-[3(log/(c C )) a! -2(logx) a ]} 
3 (logf{x) 



(3.1) 



exp < — 2(logx) c 



1 



logx 

log(2a(log x) a ~ ) 
logx 



Since the exponent in (|3.ip converges to — oo as x — > oo, Assumption [5] is satisfied and this pair 
(X,Y) satisfies the Assumptions in Section [2.1.11 

Example 3.2. Suppose, X and Y are independent and identically distributed with common dis- 
tribution F, where for a > 1, 

exp(— (logx)") if x > 1, 
1 if x < 1. 



F(x) 



As in Example 13. 1| one can check the subexponentiality condition (|2.19|) with L = 1 and by 
Corollary 12.91 F is subexponential. Hence, 

P(X + y >x) ~2P(X >x). 

Example 3.3. Suppose, X ~ Lognormal(/i, ct 2 ) and Y = e 2fl /X so that X = Y. We check the 
Assumptions in Section [2.1.11 for the pair (X, Y). The distribution Lognormal(//, <r 2 ) belongs to the 
maximal domain of attrac t ion o f the Gumbel distribution and its mean excess function e(x) has 
the form Embrechts et all 119971 . page 147, 161] 



a 2 x 



e[x 



-(l + o(l)). 



log X — (J, 

Also, (I2.ip is obvious and so, Assumption [1] of Section 12.1.11 is true. Following Remark l2.2H |3j) and 
the form of e(x), we may assume the auxiliary function 



/(*) 



a 2 x 



log x — fl 
To verify Assumption [3l fix t > 0, and note as x — > cxo, 

P(X > x, Y > tf(x)) P(X > x, e^/X > tf(x)) 







since f(x) 
x — > oo, 



P{X > x) P(X > x) 

oo. Assumption H] is verified similarly. For Assumption choose L = 1 and as 



P(X > /(x),y > /(a)) __ P{X > f(x),e 2 »/X > f(x)) 

P(X >x,Y >x) 



0. 



P(X > x) 
We conclude by Theorem 12.81 

P(X + Y > x) ~ 2P(X >x). 

Example 3.4. Example 13.31 is a special case of a more general phenomenon. Suppose, F G 
MJ9j4(A) with auxiliary function f{x) having the property 

(3.2) liminf/(x) = «5 > 0. 

Assume that the support of F is a subset of [0, oo) and xq = sup{x : F(x) < 1} = oo and also 
that xi = inf{x : F(x) > 0} = 0. Distributions satisfying these conditions include the exponential, 
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gamma, lognormal. Define X = F*~~(U), and Y = F*~(l — U), where U ~ Uniform(0, 1). This pair 
(X,Y) satisfies the Assumptions in Section [2.1.11 

Checking Assumption [2] is easy since X and Y are identically distributed. To verify Assumption 
[3l fix t > and define et = F(^). Since, x\ = 0, we have ej > 0. Then, for large x making 
f(x) > (5/2, we have 

P(X>x,Y >tf(x)) __ P{U > F{x), 1-U> F(tf{x))) 
P(X >x) ~ P(X > x) 

P{U > F{x),l-U >e t ) P(U> F{x),U <l-e t ) 



< 



P(X > x) P(X > x) 

since F(x) — > 1, and xq = oo. Assumption [J] is similarly verified. To verify Assumption [5j choose 
L such that F(k^) > | and for x sufficiently large, 

P(X > Lf(x),Y> Lf(x)) < P(X>^,Y> M) 



P(X >x) ~ P(X > x) 

P([/>F(f),l-[/>F(f)) _ 
P(X > x) 
Hence, (X, Y) satisfy the Assumptions of Section [2.1.11 and by Theorem! 

P{X + Y >x) ~2P(A >x). 
In this example, if lim x ^oo f(x) = oo, we do not need the condition x\ = 0. 
Remark 3.5. Note, in the previous two examples a comonotonic dependence structure is used. 

Example 3.6. Suppose, X = exp(Ai), Y = exp(X2), where (X\,X2) is bivariate normal with 
correlation p £ [—1,1). For simplicit y, assume each Xj has mean p and var iance a 2 > 0. This 
example is extensively considered in lAsmussen and Roias-Nandavapal [2003]. We have already 
considered the case p = — 1 in Example 13.31 so here we consider p £ (—1,1). 

Assumptions [1] and [2] of Section f2. 1.11 are easily verified. Following the same reason as in Example 
we take the auxiliary function to be 



m 



a 2 x 



log x — p 



log fix) — a io S(iogx-/x) A* logx — p 1, /logx — p 

■log 



Observe, 



a a a a ~ \ a 2 

For Assumption [3] , we have for t > 0, as x — ► oo, 

P(X >x,Y> tf(x)) P(Xi > log x, X 2 > log tf(x)) 



< 



P(X>x) P(X 1 >\ogx) 

P(X 1 + X 2 > logx + \og(tf(x))) 
P{X X >logx) 

* ( ^(l+P) (1 ° g X + I0g(t/(X)) " 2/i) 



$ 



logx— fi > 
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• logx-tj, , log /Q)-A* 



Q 1 / ioga;-/j , log j (X)-ij. iogt\ 

^r]ogx-fj,\ 



$(l2S£^i£) 



0. 



where we used (13.3|) and the fact that $ £ MDvl(A) and therefore $ is — oo-varying [Resnickl . il 98 71 . 
page 531. Note, p < 1 entails , 2 > 1. 

J r V 2 ( 1 +p) 

For Assumption [5J choose L = 1. As cc — ► oo, we have using (|3.3p . 

P(A > /(ar),y > f(x)) __ P(X 1 > logf(x),X 2 > log /(a)) 
P(X>x) P(Xi>logx) 

$ , 2(log/(a)-/ J ) 



< P(Xi + X 2 >21og/(x)) _ V v/2 ff2 (!+/>) 



P(X 1 > log a;) ^/ (logx-/^) . 

$f 2 (l°££^)( 1+o(1 )^ 



*( 



log£— £4N 



Example 3.7. Let Xi and X 2 be independent and identically distributed with the common distri- 
bution H E MDA(A), having auxiliary function /i(-) satisfying (|3.2p and infinite right end point. 
Also, suppose, F £ MDA(A) with auxiliary function /2(-)i concentrates on [0, oo) and satisfies the 
conditions in Example 13.41 Define X and Y as 

X = F*-(U)AX 1 , Y = F-(1-U) AX 2 

where U is a uniformly dist ributed r andom variable on (0,1) which is independent of (Xi,X 2 ). 

From Proposition 1.4 of Resnickl . 1 198 71 . page 43], the distribution of X belongs to the maximal 
domain of attraction of Gumbel with auxiliary function 

h(x)f 2 (x) 



m 



fi(x) + f 2 (x) 



Hence, 



and thus, 



lim sup < lim sup -r7~T + nm sup < cxd, 

x^oo J(x) x^oo fl{X) x->oo h{X) 



liminf f(x) > 0. 

x— >oo 

Also, note, 

P(X > x) =P(U > F(x), X 1 >x) = P(U > F(x))P(X 1 >x) = F{x)H{x) 

and 

P(Y > x) =P(l -U> F(x),X 2 > x) = P(l - U > F{x))P{X 2 > x) = F(x)H(x). 
Arguing as in Example \'6A\ we can show that the pair (X, Y) satisfy the assumptions in Section 

EXU 
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Example 3.8. Here is an example of a distribution for (X, Y) where our assumptions are not 
satisfied, but the asymptotic behavior is the same as in Theorem 12,81 Suppose, X and Y are iid 
with common distribution F, where 

F(x) = exp(-x a ) a €(0,1), x > 0. 

This distribution is extensively studied in iRootzenl [19861 ] and satisfies F G MDA(A) n S. Since it 



is subexponential, 

P(X + Y >x)~2P(X >x). 

However, this distribution does not satisfy Assumption [5] of Section [2. 1.11 
Since F is a Von-Mises function, we may take the auxiliary function to be 

f(x ) JM xl ' a 

n ' ~ F'(x) ~ a ■ 
Assumption [5] is not satisfied for any L > 0, since for any L > 0, as x — > oo, 

P(X > Lf(x),Y> Lf(x)) __ [F(Lf(x))] 2 _ exp(-2 [£/(*)]< 



P(X > x) F{x) exp(-x a ) 

exp(-2(^) Q ^( 1 -")) / L« _-\ 

- exp x (1 — 2(— ) x ) \ —> oo. 



exp(— x a ) \ a 

This also shows the criteria (|2.19p for F G S is sufficient but not necessary. 

4. Linear combinations of random variables with non-negative coefficients 

This section studies linear combinations of risks X, Y with non-negative coefficients. We consider 
two cases: (i) the distributions of X and Y are tail-equivalent, and (ii) the distributions of X and 
Y lack tail-equivalence. We explicitly give the asymptotic tail behavior of the linear combinations 
of risks in the tail-equivalent case and also in one special case where tail-equivalence is absent. We 
note that one cannot expect similar behavior in the two cases. 

4.1. Tail-equivalent cases. 

4.1.1. Linear combination of two random variables with non-negative coefficients. 

Theorem 4.1. Assume, (U, V) is a pair of random variables which satisfy Assumptions [7], [3J 
cmd[5| of Section \2.1.1\ Moreover, assume that Assumption^ holds in the form 

(4.1) hm P( y> x) =ce(0,oo). 

X^OO f[U > X) 

Define, S2 = a\U + a?V and ai > 0, i = 1, 2 and set m 2 = a\ V a 2 . Then, as x — >• 00, 

P{S 2 > X) ~ P(Z7 > ^-) [l{a 1= m 2 } + Cl { a 2 =m 2 }] 

We assume £/ and V are tail equivalent, i.e. the constant c cannot be and hence both the 
marginal distributions belong to MDA(A), the maximal domain of attraction of the Gumbel. If 
lirn^oo P{V > x)/P(U > x) = 0, the asymptotic behavior of P{a\U + a 2 V > x) as x — ► 00 can be 
different as illustrated in the following example. 
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Example 4.2. Assume, (U, V) are iid random variables with common distribution F, which satisfy 
Assumptions [TJ El H] and [5] of Section 12.1.11 Define the two random vectors by (U\, V\) = (U, j* V) 
and (U 2 ,V 2 ) = (U, ±V). Both pairs (U u Vi) and (U 2 ,V 2 ) satisfy Assumptions [Q EJ H and G2 of 
Section 12.1.11 For both pairs, c = 0, i.e., 

l im ^4 = and lim^4 = 0. 

x— oo P{Ui > X) x^oo P(U 2 > X) 

Since, (U, V) satisfies the Assumptions of Theorem 14.11 we have as x — > oo, 

P{Wi + 10Vi > x) = P{3U + 2V >x)~ P(3U > x) = P{3U X > x), 

and 

P(3U 2 + lOVb > x)= P(3U + 5V > x) ~ P(5F > x) = P(10y 2 > x). 
This example illustrates we cannot expect Theorem 14.11 to hold for the case c = 0. 

We now turn to the proof of Theorem 14.11 
Proof. The case a± = a 2 is resolved by Theorem 12.81 since 

P(oi(l7 + V) > x) = P{U + V > — ) ~ (1 + c)P{U > — ). 

So the interesting cases are a\ > a 2 and a\ < a 2 and for the following, assume a\ > a 2 , the other 
case being similar. 

There is nothing to prove if a 2 = 0, so assume a\ > a 2 > which makes m 2 = a\. It suffices to 
check the Assumptions in Section [2.1.11 for X = U and Y = a 2 V/a\. For this definition of X, Y, 
we have 

(4.2) lim P i F > '» = Urn ^"f"' > ' '» = Urn ^ > S^ = 0. 

V ^ x^oo P(A > z) x-^oo P(U > X) x^oo P(U > X) 

The last equ ality is t r ue fro m (|4.1[) and the fact that the tail of any distribution in MDA(A) is 
— co-varying Resnicki . 1 198 71 . page 53]. From Theorem 12.81 and (|4.2j) . we get, as x — > oo, 

P(aiU+a 2 V >x) = P(ai(U + 02V/01) > x) 

=P(17 + a 2 ^/ai > x/01) = P(X + F > x/01) 

~P(U > x/ai) = P{U > ^-) [l {ai=m2} + cl {a2=m2} ] . 

To complete the proof, the Assumptions in Section 12.1.11 must be verified for X = U and 

Y = a 2 V/a\. Assumption [TJ is assumed in the Theorem and Assumption [2] was verified in (|4.2p . For 
Assumption^ note that U G MDA(A) and suppose /(•) is the auxiliary function of the distribution 
of U. By hypothesis, for t > 0, 

(4.3) lim P(|V| >tf(x)\U>x) = 0, 

x—>oo 

and therefore, using (14.31) . 

lim P(a 2 \V\/ai > tf(x)\U > x) = lim P(\V\ > ai tf(x)/a 2 \U > x) = 0. 

x— >oo x— »oo 

Remark 12. 21 [4"j) implies we do not need to verify Assumption [H so we check Assumption [5j For 
this we have, as x — > 00, 

P(a 2 V/ai > Lf(x), U > Lf(x)) __ P(V > ai Lf(x)/a 2 , U > Lf(x)) 
P(U > x) P{U > x) 
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< P(V > Lf(x),U > Lf(x)) 
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0. 



P(U > x) 
This proves the case a\ > a<i- □ 

4.1.2. Linear Combination of more than two random variables with non-negative coefficients. 

Corollary 4.3. Assume, X\,X2, ■ ■ ■ X d are non-negative random variables which pairwise satisfy 
Assumptions [3|, E[ [31 o/ Section T2.1. 1\ Further suppose the distribution of X\ satisfies Assumption 
[7] of Section \2.1.1\ and that 

P(Xi > x) 



(4.4) 

Set c\ 

define, 



Then 



lim 



Ci G (0,oo) 



1,2, 



,d. 



x^oo P{X\ > x) 

1 and define for d > 1, S d = a\X\ + 02-^2 + • • • CLdXd, for ai > 0, 



1,2, 



, d. Also, 



m d 






and N d 



J2 <*■ 

{l<i<d:ai=m l i} 



P(S d > x) ~ iY d P(X! > 



X 



00. 



m d 



This result is co nsistent with the case wh ere X\ , X2 , • • • , X d are iid with common distribution in 
MDA(A) n5; see lDavis and Resnickl |l988^ . 

The random variables X\ , X2 , ■ ■ ■ X d are tail-equivalent and satisfy Assumption [3] of Section 
12.1.11 pai r wise. Therefore Remark 12. 2t [2|) implies pairwise asymptotic independence and hence, by 



Resnickl , Il987l . page 291], X\ . . . , X d are asymptotically independent. 



In the special case that the random variables are identically distributed, Na = \{1 < i < d : a\ = 
mrf}|, where | • | is the size of a set. 

Remark 4.4. It is possible to prove Corollary 14.31 using Corollary 12.101 However, in the proof 
it is usually difficult to verify Assumption 0] of Section 12.1.11 Note, a similar problem is avoided 
carefully in the proof of Theorem 14.11 through the help of Remark I2.2HJ41) . Though a remark similar 
to Remark I2.2IH1) could also be made for Corollary 12.101 it is notationally inconvenient. So, to 
avoid this notational difficulty, Theorem 14.11 is used for the proof. 

Proof. Proceeding by induction, note the base case for d = 2 is proved in Theorem 14.11 As an 
induction hypothesis, suppose the result is true for d = k, so, as x — ► 00, 

N k 



P(S k > x) ~ N k P(X l > — ) 

m k c k+1 

To prove the result for d = k + 1, notice, 

(4.5) mjfe+i = mfc V Ofc+i, 

and 

N k+1 = 
so that 

A^fe+i 



P(X k+1 > —) 
m k 



[ c k+l i -{a k+1 =m k+1 } + ^k i -{m k =m k+ 



1} ]> 



(4.6) 



Cfc+l 



{a k+1 =m k+1 } "I" 



N k 

Cfc+l 



{m k =m k+1 } 
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By the induction hypothesis, 



P{m~ l S k > x) __ Pim^Sk > x) P(X X > x) N k 



(4.7) lim — - — = lim 

x^oo P(X k+l > x) x^co p(Xi > x) P(X k+1 > x) c k+l 

If we prove the assumptions in Theorem 14.11 are valid with U = X k +\ and V = m^ S k , then, 
Theorem |4"7T] ([43]) . P~B]) and ([477]) imply, as x ->■ oo, 

P{S k+ i > x) = P(a k+1 X k+1 + m k S k > x) 

N k+1 



r<*j 



-P(X k+1 > ) ~ N k+1 P(X 1 > 



Cfc+i rn k+1 m k+ i 

and by induction, our result holds for all d > 2. 

Assumption [1] is assumed. For (|4.ip . consider that on the one hand, 

k 

(4.8) N k = J2 c *^ A Ci > ° 

{l<i<fc:ai=mj.} i=l 

and on the other, 

(4.9) iV fc = J] q < k \/ a < oo, 

{l<i<fc:ai=mj.} i=l 

and therefore the limit in (|4.7p satisfies N k /c k +i £ (0, oo). 

Next, suppose, two random variables [/ and V are tail equivalent and both belong to MDA(A). 
If /(•)> /(•) are the auxiliary functions of U and V respectively, then f(x) ~ /(^), as x — ► oo; 



sec 



Resnickl |1971b| ja| . Since, in the present case, all the random variables are tail-equivalent, 



Remark [22U3D implies we can work with the auxiliary function of any one of them, say X k+ \. So, 
X\ , X2 , • • • , Xd satisfy Assumptions [4] and of Section 12.1.11 pairwise with the auxiliary function 
/(•) of X k+ \. That is, for 1 < i ^ j < d, and any t > 0, 

(4.10) lim P(Xj > tf(x)\Xi >x) = 

x—>oc 

and for 1 < i < j < d, for some Lj,- > 

u in P(X i >L ij f(x),X j >L ij f(x)) = 

K ' ' P{X t > x) 

To verify Assumption El observe for t > 0, that as x — ► 00, 
PGmj^Sfcl > t/(x)|A fc+1 > x) 

< P(aiXi + a 2 X 2 H h a fc X fe > m k tf(x)\X k+1 > x) 

k k 

<Y, P ( X i > ar l m k tf(x)/k\X k+1 > x) < J2 P ( X * > tf(x)/k\X k+1 > x) -► 0, 
i=l i=l 

by (|4.10|) . For Assumption HI note, 

Pim^Sk > x) 
4-12 lim p/L * j = N k , 

x->oo P(Xi > X) 

and for 1 < i < Ar, 

,. P{{m~ l S k > x) n (m: 1 ^ > x)) P(m7 l aiXi > x) 

lim ; ^ = hm 



2^00 P(Ai > x) z^oo P(Ai > x) 



lim 
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( 4 -!3) = Cil {ai=nlk} . 

The first equality uses the assumption that X^s are non-negative. The second equality is true from 
(|4.4j) and the fact that the tail of any distribution in the maximal domain of attraction of Gumbel 
is —oo- varying. Now, for 1 < i < j < k, using (|2.5p . 

P((m7 §k > x) n (m7 aiXi > x) n {m7 cljXj > x)) 

PjXl > x) 

P((m7 ciiX; > x) n (rar o,X > x)) 

< lim - *• — — — -f — y —^ — J — J - 

z^oo P(Xi > x) 

< Um p((Xj >x)n (x, > X )) = Q 

z^oo P(X 1 > x) 

Therefore, using (|4.13p . 

Pdm^Sk >x)D (U^im^OiXi > x))) 



lim 

x^oo P(Xi > X) 

uu) ,■ Eti ^(K" 1 ^ > x) n (m^q t X t > x)) 

(414) " £& P ( Xl > X ) = Nk - 

From (14121) and (|4"l4"j) it follows that 

P((m^g fc > x) n (uJi^V^j > ^)) C ) _ n 

™ P(Zi > x) 

and this, along with (|4.4p and (|4.7p give 

(4 15) Hm P((m7 l S k > x) n (Uj^m^a,;^ > x)) c ) = q 

™ P(m7 l S k > x) 

Now, we check Assumption [U For £ > 0, as x — ► oo, 

p,,v |^,,, m -ia . x P(A" fc+1 > tf(x),m7 l S k > x) 

P(\X k+1 \ > tf(x)\m k S k >x) = — 

P(m fc Sfe > x) 

^ P(X fc+ i > tf{x),m7 1 S k > x^^jm^aiXi > x}) 

Pim^Sk > x) 

P(X k+1 \ > tf{x),dl =1 {m7 1 a i X i > x}) 



< 



< 



P(S k > m k x) 
St=i -P(^fc+i > tf(x) 1 m7 1 a i Xi > x) 
P(S fc > m fc x) 



where we have used (|4.15p . Using our induction hypothesis, we get that the quantity above is 
aymptotically equivalent to 

EiLi P ( X k+i > tfix^m^aiXi > x) 



N k P(X 1 > x) 
< EL 1 P(X k+1 >tf(x),X l >x) 



N k P(X 1 > x) 
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Eti P(Xk+i > tf(x), X t > x) P{X t > x) 

P(Xi > x) N k P(X 1 > x) ~* ' 

by gTUD . 

For assumption let, L = kL max , where L max = maxi<j<fc Li ;k +i (recall, equation (|4.1ip ) . 
Then using (|477|) . (|4TT5j) and ([4TTT]) . we have 

-P(X fe+1 > kL max f(x) 1 m^ 1 S k > kL max f{x)) 
P(X k+l > x) 
P(X fc+ i > kL max f(x),m' l : l S k > x,U , * =1 { m - 1 ai x i > x}) 

P(X fc+1 > x) 
P(X fc+ i > kL max f(x),U^ =1 {m^ 1 a i Xi > L max f(x)}) 
P(X k+1 > x) 

T,i=i p ( x k+i > L iM if(x),m~ l : 1 aiX i > L i|fc+1 /(x)) 
P{X k+1 > x) 

Eti^fc+l > L itk+1 f(x),X i > L iM1 f(x)) 

P(X k+1 >x) "* ' 



< 



< 



< 



D 



4.2. One special case where the distributions are possibly NOT tail-equivalent. 



Theorem 4.5. Assume, Yi,Yz, ■ ■ -Yd are identically distributed non-negative random variables. 

Also, assume aj, Pi > 0, i = 1, 2, . . . , d. For d > 1, define, §d = aiYf 1 + 02^2 + . . . + o-dX^ and 

set 

d 

P = \J Ph 1d= V <H, 

i=l {I<i<d-Pi=f3} 

Jd = \{l<i<d:P i =p,a i =q d }\ 
where \ ■ \ denotes the size of the set. Suppose, qdY{ , qdX 2 , ■ ■ ■ 1dY d pairwise satisfy Assumptions^ 

2] and [5| of Section \2.1.1\ and that the distribution of qdY{ satisfies Assumption^ of Section \2.1.1\ 
where the auxiliary function f(x) satisfies the additional condition that f(x) — > oo, as x — > oo. 
Then, 

P(S d >x)~ J d P( Y f > -) 

Qd 

Remark 4.6. If j3\ > 02, then Yf 1 and Y 2 are not tail-equivalent. Note, in this case, the 
asymptotic approximation of P(a{Y^ 1 + 02^2 > x) does not depend on 02- 

Theorem 14.51 shows different tail behavior from the tail-equivalent case s but follows the paradigm 

that o nly the heaviest tails matter. The Theorem shows that Theorem 1 of lAsmussen and Rojas-Nandayapa 
20081 ] is a special case of a more general phenomenon. Let (Xi, X2, . . . , X d ) ~ N(0, S), where 



S = (pij), pa = 1, V i, pij < 1 < i < j < d. 

Let, (Y 1 ,Y 2 , . . . , Y d ) ~ (exp(Xi), exp(X 2 ), . . . , exp(X d )). Clearly, 

aiYf 1 ~ Lognormal(logai,/3f) 

From Example 13.61 {idXx 1 QdYj f, ■ ■ ■ , QdY^ ) satisfy the Assumptions of Theorem 14.51 where q d , (3 
have the same meaning as in Theorem 14.51 Also, [Z\, Z2, ■ ■ ■ , Z&) = (a\Y{ , a 2 Y 2 , . . . , a n Y^ d ) 
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satisfies the Assumptions of Theorem 1 of lAsmussen and Roias-Nandavapal 20081 ] . The results of 
that theorem and Theorem 14,51 match. 



Proof. Without loss of generality, assume (3\ = (3 and a\ = q d . Also, assume en > for i = 
1,2, ... ,d. Denote, 

X i = a i Yf i i = l,2,...,d. 

To start, suppose, for some i G {2, . . . , d}, fy < j3. Then, for large x, [aiYf 1 > x] C [^fYf > x], 
and hence for large x, 

P(aiYf* > x) < P(^Yf >x) = P(q d Yf > 2x). 

Then, 

(4. 16) c, = li,„ gf > ') = , im EMtll> < lim Wf>M _ „ 

x-*oo P{X X > x) x-,00 pfaYf > x) x ^°° P(q d Y[ > x) 

Next, suppose, for some i € {2, . . . , d}, fy = (3, en < q d . Then, 

(4.17, c, = li m %*>*> = lin, P <°-^ > " = li m ^ > gj = 0. 

x^oo P(Xi > re) x^oo p( qd Y{ > x) *-°° Pfe^T > x) 

In both the equations (J4.16P and (|4.17p . the last equalities are true from the fact that the tail of 
any distribution in the maximal domain of attraction of the Gumbel is — oo- varying. 
Finally, suppose, for some i € {2, . . . , d}, /3j = f3, Oj = qd- 

, A ,K v g(gj > x) P(Yf > f d ) P(Yf > f d ) 

(4.18) Cj = lrm —7— = hm g — = hm 5 — = 1. 

z->oo P(Xi > x) a;-*oo ply? > — ) x - + °° PfY^ > — ) 

It suffices to check the assumptions in Corollary 12. 101 with this set of X\, X2, ■ ■ ■ , X d , since then 
Corollary ElO] and (l4~16l) . (liTfl) . (|4~T8l) would imply, as x -> 00, 

d 
P(S d > x) ~ (1 + £ q)P(A x > x) ~ J d P(Xi > x) = ^^(y/ 3 5 

i=2 

Assumption [1] is assumed in the Theorem statement and (|2.22p is already shown in (|4.16p . (J4.17J) 
and (|4.18p . For assumptions [3] and HI proceed as follows. By hypothesis, we know that X\ belongs 
to the maximal domain of attraction of the Gumbel distribution. Let /(•) be the auxiliary function 
corresponding to the distribution of X\. By hypothesis, we know, for t > 0, for 1 < i 7^ j < d, 

(4.19) lim P(q d Yf > tf(x)\q d Y? > x) = 0. 



x_ s 

q d ' 



x— »oo 



Using Remark I2.11l |3|) , it is enough to show 

(4.20) lin F(X J -> t /(x),X i >x) = 

x^oo P\X\ > X) 

and to see this, note that since /(x) — ► 00, for large x and for all t > 0, [X,- > tf(x),Xi > x] C 
[ gd F/>t/(x),^y/>x]. Hence, 

f >tf(x),q d Yf 
P{Xi > x) " *££> P(g d l? > x) 



,. P(X,- > t/(x), X > x) P(g d 3f > tf(x),q d Yf > x) 

hm — t— : < hm - 3 = 0. 



x— >oo 
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For Assumption [5j using Remark 12. lll |3|). we show, for some L > 0, 

P(X 3 >Lf(x),X i >Lf(x)) 



lim 

x— >oo 



0. 



(4 ' 21> -•--'- P« > x) 

By hypothesis, we know, for all 1 < i < j < d, there exists some Ly > 0, 



(4.22) 



lim 

x— »oo 



P( gd y/ > Lijfix^qgYf > Ljjfjx)) 

P(q d Yf > x) 



0. 



Also, note that since /(jc) -> oo, for large x, [Xj > Lijf{x),Xi > Lijf(x)] C [%*/ > Lijf(x),q d Yf > 
Lijf(x)]. Hence, 



lim P(Xj > bilMl^l > hilMl < l im PJQdYf > Lij/ixlggYf > Ly/(x)) 



P(Ai > x) 



P(^F 1 P > x) 



by g22D. 



D 



5. An Optimization Problem 



5.1. The problem. Suppose, we have a portfolio consisting of d financial instruments. The risk 
per unit of the i-th instrument is JQ. The goal is to earn revenue $L. Assume, each unit of the 
z-th instrument earns $/j over the chosen time horizon. Subject to earnings being at least $L, how 
many units of each instrument, 01,02, .. . , «d, should be used to build the portfolio, so that the 
probability that the total portfolio risk 01X1 + 02X2 + . . . -\-adXd exceeds some fixed large threshold 
x, is minimal? 

Thus, consider the following optimization problem: 

d 

y^ ai Xi > x 

s.t. aili + a 2 /2 H 1- add > L, 

Oj > 0, i = 1, 2, . . . ,d. 
For a more general case, consider the following optimization problem: 

" d 

y diXi > x 



min P 

{a!,...,a d } 



min P 
{ai,...,a d } 



i=\ 



s.t. h(ai,a 2 ,... ,a d ) > L, 
a.j > 0, i = 1, 2, . . . , d. 

5.2. The method. Suppose, Xi,X2, . . . ,X d satisfy the assumptions of Corollary 14.31 Even with 
these assumptions, exact solution of the optimization problem is difficult. An obvious way to obtain 
an approximate solution to the optimization problem is to assume that the threshold x is big and 
use the asymptotic approximation of P(a±Xi + 02X2 + . . . + a d X d > x) from Corollary 14.31 hoping 
that the solution of the resulting optimization problem is close to the actual optimal value. So, 
using the notation of Corollary 14.31 we solve the following optimization problem: 

min N d P(X 1 > JL) 

{ai,...,a d \ m d 

s.t. /i(ai,a 2 ,... ,a d ) > L, 
a% > 0, i = 1,2, . . . ,d. 
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Suppose ai, a,2, ■ ■ ■ , o d and di,a,2, ■■ ■ ,a d are two feasible solutions for the given set of constraints. 
Set 

d 
m d =\f Oi, N d = ^2 a 

i=l {l<i<d:a,i=m d } 

d 
m d =\J a,i, N d = Y^ °i 

i=l {l<i<d:ai=m l i} 

If rhd > m d , then since, P[X\ < x] E MDA(A), as x — > 00, 

P(X! > x/m d ) 

— ^ CXD 

P(Xi > x/m d ) 
Now, since both TV^, iVrf S [Af =1 Cj, d Vf =1 Cj], we have as x — > 00, 

iV d P(Xi > x/rhd) 



N d P{X 1 > x/m d ) 



00. 



So, we hope that 0,1,0,2, ■ ■ ■ , a d is a better feasible solution for the optimization problem. 

Thus, values of 01, 02, ■ ■ ■ , a d which solve the above optimization problem can be computed by 
solving the following two optimization problems in sequence. 

(i) First solve 

min m d = max{ai, 02, ■ ■ ■ , a d } 

{ai...,a d } 

s.t. h(ai,a 2 , ... ,a d ) > L, 
a>i > 0, i = 1, 2, . . . ,d. 

(ii) Suppose, the best choice of a±, 02, ■ ■ ■ , a d gives m as the value of the objective function for 
the optimization problem in (i). Then we solve 

min N d = V, c i 

{ai,...,a d } r . ..fr , 

{±<i<d:ai=my 

s.t. h(ai,a 2 ,... ,a d ) > L, 

max{ai,a2, • • • ,a d } = m, 
ai > i = 1,2, ... ,d. 

5.3. A special case. The motivating case is that h is a linear function with positive coefficients 
of the form 

h(ai,a 2 ,...,a d ) = a\h + a 2 h + ••• + a d l d . 
The approximate solution using the asymptotic form of P[%2i=i a iXi > x] is 

a 1 = a 2 = ■■■ = a d = L/(h + h + ■■■ + Id)- 
This leads to m = L/{1\ + I2 + ■ . . + Id) and N d = ^2 i= i ci. 
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6. Simulation studies 

We carried out some simulation studies to check for fixed large thresholds the accuracy of the 
asymptotic approximation in Theorem 12.81 and also to check how good is the approximate solution 
for the optimization problem. As expected, in some cases the approximation works well whereas in 
others it performs poorly which suggests caution about using the asymptotic results for numerical 
purposes. Simulation also suggests that the approximate solution of the optimization problem 
works well in cases where simulation studies suggest that the approximation is good for fixed large 
thresholds. One particular model studied, Example l3.6l with fj, = 0, a = 1, is noted here to illustrate 
the point. We varied p and observed the asymptotic behavior of the sum of the risks. 

6.1. Where is the approximation good? To test the approximation for P(X + Y > x), we have 
to find good simulation estimates of the probabilities P(X + Y > x). This, however is not easy, 
especially in the case when the marginal distributions of the risks X and Y are sub exponential 
and is still a topic of current research in the simulation community. The approach usually taken 
in these cases is Conditional Monte Carlo Asmussen and Glvnnl . 120071 . page 173]. So, this method 



is used to compute P(X + Y > x) and the simulation estimates are compared with the theoretical 
approximations. 



The simulation of P(X+Y > x) uses the algorithm suggested in l Asmussen and Roias-Nandavapa 



20081 ] for p G (—1, 1) who also note the properties of this algorithm. If p = —1, we have a way to 



compute the probability exactly. In this case, X = 1/Y almost surely, so in the following manner 
we compute the required probability: 



PIX + 1>X)=P(X> X + V ^ ) + P(X< X -^ 



X J \ 2 / V 2 



- + ^ g ~* ))+p(logX<log / 



x + v ' x 2 — 2\ \ / / x — \fx 1 ~ 



$ log +$ log 



6.1.1. Patterns in the results. For judging the quality of the asymptotic approximation, we focus 
on the simulation estimate P(X + Y > x) and not the threshold x, since a change of distribution 
may imply a change in how rare is a particular threshold crossing. So, when comparing the quality 
of the asymptotic approximation across different models, it makes more sense to focus on the value 
of P(X + Y > x), rather than the particular threshold x. When p = —1, exact calculations suggest 
that the approximation is extremely good even when the actual probability P(X + Y > x) is of 
the order of 10~ 2 . As expected, the asymptotic approximation improves as a function of increasing 
threshold. When p S (—1, 1), we rely on the simulation estimate as a surrogate for the exact tail 
probability and compare it with the theoretical approximations. 

The results indicate that the closer p is to —1, the better the approximation. For p = —1, the 
approximation is good for events with probability of the order of 10 -2 and to achieve comparable 
precision in the relative error when p = 0, the event has to be much rarer and have a probability of 
the order of 10 -10 . For p = 0.9, the results for different thresholds did not show any convergence 
pattern. This emphasizes that in practice the numerical approximations should be used with 
caution. Clearly for p = 1 the asymptotic approximation is not correct and p = 0.9 is expected to 
behave somewhat like the case when p = 1. 
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Table 1. p = — 1 
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Threshold 


Actual probability 


Asymptotic approximation 


Ratio 


10 


0.0219 


0.0213 


1.0272 


16 


0.0056 


0.0056 


1.0121 


24 


0.0015 


0.0015 


1.0060 


30 


6.7365 xl0~ 4 


6.7091 xl0~ 4 


1.0041 


100 


4.1233 xl0~ b 


4.1213 xl0" b 


1.0005 


1000 


4.9238 xlQ- vz 


4.9238 xl0~^ 


1.0000 



The tables give representative results. We first give the results for p = —1 in Table [U since in 
this case no simulation is required. The column 'Ratio' in Table Q] is defined as 

Actual probability 

Ratio = — . 

Asymptotic approximation 

For subsequent tables, the columns 'Ratio' and 'Half-width' are defined as 

Simulation estimated probability 

Ratio = 

Asymptotic approximation 

Half-width =Half-width of the 95% confidence interval of the ratio. 
In each case, 10 observations were used to compute the probability estimates. 



Table 2. p 



-0.9 



Threshold 


Simulation estimated probability 


Asymptotic approximation 


Ratio 


Half-width 


3 


0.3687 


0.2719 


1.3556 


0.0006 


5 


0.1207 


0.1075 


1.1227 


0.0012 


10 


0.0221 


0.0213 


1.0375 


0.0026 


20 


0.0028 


0.0027 


1.0082 


0.0064 


30 


6.8873 xl0~ 4 


6.7091 xl0~ 4 


1.0265 


0.0119 


40 


2.2134 xl0~ 4 


2.2524 xl0~ 4 


0.9827 


0.0183 


50 


9.3675 xl0~ b 


9.1526 xl0~ b 


1.0235 


0.0285 


Table 3. p = 


Threshold 


Simulation estimated probability 


Asymptotic approximation 


Ratio 


Half-width 


10 


0.0338 


0.0213 


1.5844 


0.0033 


50 


1.0798 xl0~ 4 


9.1526 xl0~ 5 


1.1798 


0.0002 


100 


4.5032 xl0" b 


4.1213 xl0" b 


1.0927 


0.0001 


300 


1.2117 xlO" 8 


1.1718 xlO" 8 


1.0341 


0.0000 


600 


1.6147 xl0~ iu 


1.5853 xl0~ iu 


1.0185 


0.0122 


1000 


4.9821 xl0~ 12 


4.9238 xl0~^ 


1.0118 


0.0000 


2000 


1.9620 xlO -14 


2.9310 xl0~ i4 


1.0106 


0.0000 



26 



A. MITRA AND S. I. RESNICK 

Table 4. p = 0.9 



Threshold 


Simulation estimated probability 


Asymptotic approximation 


Ratio 


Half-width 


10 


0.0521 


0.0213 


2.4439 


0.0088 


30 


0.0030 


6.7091 xl0~ 4 


4.4081 


0.0275 


50 


5.2652 xlO" 4 


9.1526 xl0~ 5 


5.7527 


0.0759 


75 


1.1217 xlO" 4 


1.5781 xl0" b 


7.1077 


0.1843 


100 


3.4333 xlO - " 


4.1213 xl0" b 


8.3307 


0.3642 



6.2. How good is the portfolio suggestion? Here, we consider the quality of our approximate 
solutions for the optimization problem. We choose the same risk model given in Example 13.61 
because we have information about which values of p lead to good asymptotic approximation. We 
resort to a naive method for analyzing the optimization. For different (ai, 02), we obtain estimates 
of P(a\X + CL2Y > x) through simulation. To get the estimates proceed as follows: For 01, a 2 > 



a ± X 
a 2 Y 



exp{log(ai) + Xi} 
exp{log(a 2 ) + X 2 } 



pe[-M) 



Now, 

' Z 1 \ __ ( \og{a 1 )+X 1 \ (( log(a x ) W 1 / 

Z 2 J~\ log(a 2 )+X 2 J \\ log(q 2 ) )'\p 1 

So, again we are in the framework of lAsmussen and Roias-Nandavapal 20081 ]. and we use the 
algorithm given in their paper to estimate the rare event probabilities. When either a\ or a 2 
is equal to 0, we can compute the exact probability and hence do not need an estimate. We 
choose (0,1,0,2) in the following way. Let C be the set of all possible (01,02) which satisfy the 
constraint. First, a\ is chosen from the corresponding projection of C with a small grid, and then 
for each a±, a 2 is determined from the constraint. Let us call this set C*. For (01,02) G C*, 
P(a\X + 02^ > x) is estimated through simulation and then it is observed which (01,02) gives 
the minimum estimate of P(a±X + a 2 Y > x). Let, (5i, a 2 ) be this pair; i.e. P(a\X + a 2 Y > x) = 
min( ai>a2 ) G (7» P(a\X + a 2 Y > x). Also, let (a*,^) be the approximate solution of the optimization 
problem as noted in the previous section. Relative error of the approximate solution is computed 
by comparing P(a\X + a 2 Y > x) with min( aia2 ) eC » P(a\X + 02^ > x). 

6.2.1. Identifying patterns. We do not have error estimates for our simulation results. One could 
consider bootstrapping to obtain such error estimates, but we have not done so. Despite the 
weaknesses of the naive procedure, the results are interesting. 

We note one case with the linear constraint 2ai+3a2 = 1. The suggested optimum portfolio based 
on asymptotic approximation is (o*, a|) = (0.2, 0.2). The 3 cases where p = —0.9, 0, 0.9, are chosen, 
the reason being that we know from the results in earlier section that the asymptotic approximation 
is good in the case p = —0.9, reasonable when p = and rather bad when p = 0.9. The approximate 
solution (a*,a 2 ) relies on replacing the original objective function by its asymptotic approximation, 
and so it is reasonable to expect different accuracies for these three values of p and this turned 
out to be the case. In the cases of p = —0.9 and p = 0, we see that 5i comes close to 0.2 as the 
threshold x increases. But, in the case of p = 0.9, no pattern in the convergence of a\ is observed 
which is expected because for p = 1, both the risks are actually the same random variable which 
implies indifference to the choice of (01,02) E C. 

Another remark is that in each case of p = —0.9, 0, 0.9, the relative errors do not show any 
convergence pattern. Perhaps to expect otherwise is unrealistic as we are using the minimum of 
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some simualtion estimates to compute the relative error. Still, we illustrate through an example 
the accuracy by comparing with an extreme case where we build the portfolio consisting of only one 
asset. For p = 0, and threshold x = 10, the extreme cases will yield probabilities 0.2441 and 0.1360. 
These risk probabilities are quite high compared that of our suggested optimal portfolio {a\,a^) 
based on asymptotic approximation, which has risk probability P{a\X + a^Y > x) = 1.0793 x 10 -4 ; 
also, the minimum of the simulation estimates P{a\X + a^Y > x) is of the same order. So, the 
suggested portfolio (a^a^) is quite effective in reducing the risk and possibly close to the best one. 

The following additional conclusion can be made. In the case of p = —0.9, even when P{a\X + 
CL2Y > x) is as big as 0.11, it is quite close to P(a\X + a^Y > x), indicating that the suggested 
optimal choice (aj, a^) significantly reduces the risk probability. For p = 0, a comparable statement 
can be made when the minimum of the probability estimates is of the order of 10 -2 . However, for 
p = 0.9, the relative errors are never close to 0. Interestingly, even for p = 0.9, P(a\X + a 2 Y > x) 
and P{a\X + a^Y > x) are almost always of the same order. However, it should be noted at this 
point that even in this case of p = 0.9, the extreme cases where the portfolio is built on entirely 
one of the assets, P{a\X + a 2 Y > x) is of a much bigger order than P{a\X + a 2 Y > x). So, 
in this case, possibly P(a\X + a 2 Y > x) differs considerably from choices where 01,02 > and 
the case where either a\ = or 02 = 0, but does not differ too much among the choices where 
(01,02) £ C, 01,02 > 0. This fact justifies the intuition as mentioned before that the case p = 0.9 
is similar to case p = 1 . Some of the results are noted in tables below. 

Results are summarized in the tables for p = —0.9, 0, 0.9 and constraint 2a\ + 3a2 = 1. For each 
fixed p, we give 

• the threshold x, 

• 01, where (01,02) £ C* and 



P(aiX + 02^ > x) = min P(a\X + a 2 Y > x 



El = min (ai>a2)GC .* P(a 1 X 
E2 = P(a\X + a%Y > x), 
the 'Relative error' 



02^ > x), 



E2-E1 



El 



For each value of p, ai is chosen with gap 0.01 from the projection of C*; i.e. we considered 
(ai = 0,0.01,0.02, . . ., 0.5). For each such a\, we used 10000 observations to obtain the estimates 
of the probability P(a\X + a 2 Y > x). 



Table 5. p 



-0.9 



Threshold 


Oi 


El 


E2 


Relative error 


1 


0.13 


0.1097 


0.1204 


0.0975 


3 


0.18 


0.0067 


0.0069 


0.0322 


5 


0.19 


0.0013 


0.0013 


0.0294 


10 


0.19 


1.0299 xl0~ 4 


1.0592 xl0~ 4 


0.0284 


20 


0.21 


2.0806 xl0~ b 


2.0806 xl0~ b 


1.2213 xl0~ 1& 
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Table 6. p = 



Threshold 


a\ 


El 


E2 


Relative error 


1 


0.03 


0.1349 


0.1723 


0.2765 


3 


0.16 


0.0093 


0.0101 


0.0759 


5 


0.18 


0.0016 


0.0017 


0.0503 


10 


0.19 


1.0424 xl0~ 4 


1.0793 xl0~ 4 


0.0354 


20 


0.20 


4.3888 xl0~ b 


4.3888 xl0~ b 





Table 7. p = 0.9 


Threshold 


«i 


El 


E2 


Relative error 


1 


0.01 


0.1360 


0.1798 


0.3223 


3 


0.01 


0.0140 


0.0208 


0.4831 


5 


0.02 


0.0033 


0.0050 


0.5146 


10 


0.02 


2.8357 xl0~ 4 


4.9475 xl0~ 4 


0.7447 


20 


0.04 


1.3241 xl0~ b 


2.4023 xl0~ b 


0.8142 



7. Concluding Remarks 

An important case for the study of asymptotic behavior of the sum of risks is the case where the 
risks are asymptotically independent, identically distributed and belong to the maximal domain 
of attraction of the Gumbel distribution. Many commonly occuring risk distributions fall in this 
category. We have provided sufficient conditions for 

P(X + Y>x) 



lim 

z->oc P(X > X) 

and extended the conditions to cover the case where the marginal distributions are not the same 
and to the case where some risk distributions have lighter tail but the distribution does not belong 
to the maximal domain of attraction of the Gumbel. We are not able to provide necessary and 
sufficient conditions for this kind of asymptotic behavior which is an unresolved problem. It will 
be interesting to see if it is possible to find a distribution of risks (X, Y) for which the risks 
are asymptotically independent, identically distributed, belong to MDA(A), and the asymptotic 
behavior of the sum is different than two cases mentioned in the introduction, viz. 

.. P(X + Y>x) ro . 
hm V G {2, oo} 

rr^oo F{X > X) 

Even for cases where the asymptotic behavior is understood, nothing is known about the rate of 
convergence in these cases; i.e. a quantitative estimate how good the approximation 2P(X > x) 
is for the quantity P(X + Y > x) for a large threshold x. Simulation studies indicate in certain 
circumstances the approximation is accurate, but in other cases its accuracy is dismal. 

We observed in the previous section that when tail probability approximation is good, the approx- 
imate solution of the optimization problem is also accurate whereas in the other cases this solution 
has poor accuracy. So, results on the rate of convergence would contribute to understanding the 
appropriateness of the approximate solutions in different scenarios. 

An anonymous and conscientious referee provided many insightful and helpful comments. 
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